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PREFACE 


This  investigation  was  conducted  as  part  of  the  Flood  Control  Hydraulics 
Research  and  Development  Program,  sponsored  by  the  Office,  Chief  of  Engineers 
(OCE),  US  Army,  and  administered  by  the  Hydraulics  Laboratory  of  the  US  Army 
Engineer  Waterways  Experiment  Station  (WES). 

The  study  was  accomplished  under  the  direction  of  Messrs.  H.  B.  Simmons 
and  F.  A.  Herrmann,  Jr.,  former  and  present  Chiefs  of  the  Hydraulics  Labora¬ 
tory;  J.  L.  Grace,  Jr.,  Chief  of  the  Hydraulic  Structures  Division;  and  J.  P. 
Holland,  Chief  of  the  Reservoir  Water  Quality  Branch  (Physical).  Mr.  M.  B. 
Boyd,  Chief  of  the  Hydraulic  Analysis  Division,  was  the  laboratory  Program 
Manager,  and  Mr.  Tom  Munsey  of  OCE  was  Technical  Monitor.  The  study  was 
conducted  and  this  report  was  prepared  by  Dr.  Robert  S.  Bernard.  This  report 
was  edited  by  Mrs.  Beth  F.  Vavra,  Publications  and  Graphic  Arts  Division. 

Director  of  WES  was  COL  Allen  F.  Grum,  USA.  Technical  Director  was 
Dr.  Robert  W.  Whalin. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  metric 
(SI)  units  as  follows: 


_ Multiply _ 

cubic  feet  per  second 
feet 

feet  of  water 


_ By _ 

0.02831685 
0.30M8 
2989.0  . 


_ To  Obtain 

cubic  metres  per  second 

metres 

pascals 
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PRESSURE  CALCULATION  FOR  TWO-DIMENSIONAL 
FLOW  INSIDE  HYDRAULIC  STRUCTURES 


PART  I:  INTRODUCTION 


1.  ''-When  the  static  pressure  inside  a  hydraulic  structure  falls  below 
some  critical  value,  voids  containing  water  vapor  form  therein.  The  voids 
then  collapse  with  great  force,  pitting  and  eroding  the  internal  surfaces  of 
the  structure.  The  formation/collapse  of  these  voids  is  called  cavitation, 

•/ 

and  the  prevention  thereof  is  a  prime  concern  in  structural  design. 

2.  The  prediction  of  cavitation  is  generally  accomplished  by  empirical 
means,  using  formulas  that  require  knowledge  of  the  local  static  pressure 
(Robertson  1965).  The  latter  is  obtainable  from  simple  empirical  equations  in 
conduits;  but  without  site-specific  physical  modeling,  such  equations  have  not 
been  available  for  structures  with  complex  internal  geometry.  Thus,  for  hy¬ 
draulic  structures  in  general,  one  must  solve  the  governing  equations  (the 
Navier-Stokes  equations)  in  order  to  predict  the  pressure.  Given  a  means  of 
computing  pressure  as  a  function  of  internal  geometry,  pool  elevation,  and 
flow  rate,  one  can  then  screen  ideas  for  new  configurations  and  structural 
modifications  with  regard  to  cavitation. 


3.  Discussed  briefly  in  Appendix  A  is  a  method  that  has  been  developed 
for  calculating  velocities  and  pressures  inside  a  two-dimensional  (2-D)  struc¬ 
ture  of  arbitrary  shape.  The  procedure  requires  the  numerical  solution  of  the 
Navier-Stokes  equations  (in  stream-function/vorticity  form),  in  order  to  de¬ 
termine  the  velocity  distribution.  The  pressure  is  then  computed  from  the 
known  velocity  and  vorticity  field  by  numerical  integration  of  the  momentum 
equation.  This  procedure  is  implemented  in  the  VORTEX  computer  code,  which  is 
a  modified  version  of  the  WESSEL  finite  difference  code  (Thompson  and  Bernard 
1985).  Like  WESSEL,  the  VORTEX  code  makes  hydrodynamic  computations  using 
boundary-fitted  finite  difference  grids  generated  by  WESCOR  code  (Thompson 
1983).  A  brief  discussion  of  the  theory  for  VORTEX  is  given  in  Appendix  A. 


4.  -The  VORTEX  code  has  been  used  to  calculate  pressure  distributions 
for  three  distinct  flow  conditions  in  the  outlet  works  for  Taylorsville  Lake, 
Salt  River,  Kentuefejf.  N  Results  of  these  calculations  (given  in  PART  II)  ex¬ 
hibit  favorable  agreement  with  piezometric  data  from  physical  model  tests 
conducted  by  Dortch  (1975).  Indications  are  that  the  VORTEX  code  may  be  of 
considerable  value  in  helping  to  identify  and  rectify  flow  conditions  that 
promote  cavitation. 
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PART  II:  CALCULATIONS  FOR  TAYLORSVILLE  OUTLET  WORKS 


Prototype  Description 


5.  Dortch  (1975)  reports  an  exhaustive  physical  model  study  of  the 
outlet  works  for  Taylorsville  Lake,  Salt  River,  Kentucky.  He  describes  the 
prototype  as  follows: 

The  plan  for  the  project  consists  of  a  rock-filled  dam, 
an  open  cut  uncontrolled  spillway  in  the  right  abutment, 
and  a  controlled  outlet  works  through  the  right  abut¬ 
ment.  The  top  of  the  dam  will  be  at  el  622.0*  with  the 
spillway  crest  at  el  592.0. 

Reservoir  releases  will  be  regulated  by  a  gated  intake 
tower,  consisting  of  two  flood-control  intakes  at  the 
base  of  the  structure  (el  474.0)  and  two  wet  wells  with 
five  6-  by  6-ft**  water-quality  intakes  in  each  wet  well 
at  elevations  ranging  from  503.0  to  534.0.  Both  flood- 
control  and  water-quality  flows  pass  through  two  sepa¬ 
rate  5.5-  by  14.75-ft  rectangular  gate  passages.  The 
two  gate  passages  transition  into  a  single  11.5-  by 
1 4 .75-ft  oblong  conduit.  The  last  20  ft  of  the  oblong 
conduit  contains  a  transition  to  a  flat  bottom  conduit 
before  discharging  into  an  outlet  transition  and 
stilling  basin.  A  profile  depicting  the  general  plan 
and  original  design  of  the  outlet  works  is  shown  in 
Plate  1  [Figure  1 ] . 

During  selective  withdrawal  operation,  the  emergency 
gates  will  be  closed  and  flow  will  be  discharged  through 
the  multilevel  intakes  into  the  wet  wells  and  through  an 
opening  located  in  the  roof  of  the  gate  passages  between 
the  emergency  and  service  gates.  The  service  gates  will 
be  used  to  regulate  the  selective  withdrawal  releases. 

The  locations  of  the  ten  multilevel  intakes  (five  in¬ 
takes  in  each  wet  well)  are  shown  in  Plate  2  [Figure  2]. 

An  1 8-in. -diameter  pipe  bypass  around  each  service  gate 
will  be  provided  to  regulate  the  release  of  low  flows 
with  the  service  gates  closed. 


Computer  Calculations 


6.  The  VORTEX  code  has  been  used  to  compute  the  steady-state  flow 


*  All  elevations  (el)  cited  herein  are  in  feet  referred  to  the  National 
Geodetic  Vertical  Datum  (NGVD). 

**  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 


through  the  right  half  (facing  downstream)  of  the  outlet  works.  The  flow 
region  of  interest  is  section  A-A  in  Figure  3.  Finite  difference  grids  were 
generated  for  (a)  flow  through  the  flood-control  facility  with  the  emergency 
and  service  gates  open  100  percent  (Figure  4),  and  (b)  flow  through  the  water 
quality  facility  with  the  emergency  gate  closed  (Figure  5).  In  each  case, 
only  the  anticipated  2-D  portion  of  the  flow  field  was  simulated,  and  three- 
dimensional  effects  in  the  approach  flow  and  the  upper  part  of  the  intake 
structure  were  neglected.  VORTEX  does  not  compute  flow  rate  (discharge)  as  a 
function  of  pool  elevation,  so  both  quantities  were  specified  in  advance  using 
information  given  by  Dortch  (1975).  Flow  distribution  was  assumed  uniform 
along  designated  inflow  and  outflow  boundary  segments. 

7.  For  calculating  pressure,  a  reference  point  was  chosen  on  one  of  the 
inlets,  at  which  Bernoulli's  equation  was  used  to  calculate  the  reference 
pressure: 


Pref 


pg(v 


pool 


yref 


2 

ref 


where 

yp0oi  =  P°°l  elevation 
yref  =*  reference-point  elevation 
vref  *  reference-point  flow  velocity 

The  velocity  and  pressure  at  the  water  surface  were  assumed  to  be  zero.  With 
the  reference  pressure  known,  the  differential  momentum  equation  ( A1 3 )  was 
integrated  along  a  path  lying  just  inside  the  perimeter  of  the  flow  field  to 
obtain  pressure  values  along  the  perimeter  itself  (by  extrapolation).  The 
perimeter  segment  of  greatest  interest  was  the  circular  bend  in  the  wet  well 
(piezometers  31 _35),  since  this  is  where  cavitation  was  thought  to  be  most 
likely.  The  finite  difference  grids  (Figures  4  and  5)  were  therefore  designed 
to  provide  the  greatest  resolution  along  this  bend. 


Discussion  of  Results 


8.  Figures  6  and  7  show  computed  velocity  vectors  for  a  flood-control 
operation  with  service  and  emergency  gates  both  open  100  percent.  Figure  8 
depicts  streamlines  for  the  same  hydraulic  conditions.  In  this  case  the  wet 
well  develops  into  a  region  of  recirculation,  while  the  flow  elsewhere  is 
nearly  uniform  (except  along  walls  where  the  no-slip  condition  applies). 
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Comparisons  of  VORTEX  pressure  predictions  with  piezometric  data  from  the 
physical  model  (Dortch  1975)  are  presented  in  Figures  9-11.  The  magnitudes 
and  trends  of  the  calculations  are  substantially  the  same  as  the  piezometric 
data,  though  some  of  the  details  are  lacking.  As  uncalibrated  predictions, 
however,  the  code  results  are  quite  acceptable. 

9.  Figures  12  and  13  show  computed  velocity  vectors  and  streamlines, 
respectively,  for  a  water  quality  operation  with  the  emergency  gate  closed  and 
the  service  gate  open  25  percent.  Eighty  percent  of  the  flow  enters  through 
the  top  inlet,  and  the  rest  enters  from  the  left  (Figures  5  and  12).  Regions 
of  recirculation  occur  downstream  of  all  corners,  causing  the  streamlines  to 
form  a  nearly  straight  channel  between  the  inlets  and  outlet  (Figure  13), 
P-essure  predictions  are  compared  with  piezometer  data  in  Figures  14-16.  The 
same  comments  apply  to  the  accuracy  of  these  results  as  to  those  for  the 
flood-control  flow  calculation. 

10.  A  series  of  flow  calculations  was  performed  for  operation  of  the 
water  quality  facility  alone,  with  the  emergency  gate  closed  and  the  service 
gate  open  50  percent.  The  tabulation  below  gives  the  flow  rates  and  asso¬ 
ciated  pool  elevations,  and  Figures  17  and  18  show  the  velocity  vectors  and 
streamlines,  respectively,  at  a  flow  rate  of  1,520  cfs.  Data  in  the  tabula¬ 
tion  were  taken  from  a  discharge  curve  given  by  Dortch  (1975).  Figure  19 
presents  a  comparison  of  predicted  and  measured  pressures  for  piezometer  33, 
which  registered  the  lowest  measured  pressure  along  the  circular  bend  (piezom¬ 
eters  31-35)  in  the  wet  well.  The  lowest  computed  pressure  occurred  about 
halfway  between  piezometers  32  and  33;  nonetheless,  the  agreement  between 
calculation  and  experiment  is  quite  good  for  piezometer  33.  No  other  piezo¬ 
metric  data  were  reported  for  this  series  of  tests. 

Pool  Elevations  and  Flow  Rates  with  Emergency 
Gate  Closed  and  Service  Gate  Open  50  Percent 


Pool  El 

Flow  Rate,  cfs 

530 

1,520 

540 

1,700 

550 

1 ,840 

560 

1  ,960 

570 

2,080 

580 

2,200 

590 

2,300 

PART  III:  DISCUSSION  AND  CONCLUSIONS 


11.  A  method  has  been  developed  for  calculating  pressures  for  two- 
dimensional  flow  inside  hydraulic  structures.  This  method  has  been  incorpo¬ 
rated  in  the  VORTEX  computer  code,  which  solves  the  Navier-Stokes  equations  in 
stream-function/vorticity  form.  The  code  has  been  used  to  calculate  pressure 
distributions  for  three  different  hydraulic  conditions  in  the  Taylorsville 
outlet  works.  Agreement  between  code  predictions  and  physical-model  data  was 
fair  to  good,  even  though  the  real  flow  was  probably  turbulent  and  somewhat 
three-dimensional.  No  effort  was  made  to  calibrate  the  code,  and  the  accuracy 
of  the  computed  results  illustrates  the  confidence  with  which  predictions  can 
be  made  for  the  hydraulic  conditions  simulated. 

12.  The  VORTEX  code  is  relatively  easy  to  use  after  a  few  months  ex¬ 
perience.  The  main  prerequisite  is  that  the  user  have  a  basic  understanding 
of  fluid  mechanics  in  order  to  interpret  the  results.  Previous  experience  in 
numerical  fluid  mechanics  is  helpful  but  not  essential,  since  guidance  can  be 
given  for  grid  generation  and  for  the  selection  of  code  parameters.  (The 
latter  task  concerns  numerical  stability  rather  than  predictive  accuracy.) 

The  VORTEX  code  can  be  used  by  engineers  to  screen  ideas  for  new  designs  or 
structural  modifications  at  a  modest  expenditure  of  computer  funds.  The 
steady-state  calculations  reported  herein  required  about  4,000  central  pro¬ 
cessing  seconds  on  the  CYBER  760. 

13-  No  effort  was  made  herein  to  model  turbulence  or  to  account  for 
three-dimensional  effects.  Moreover,  the  numerical  scheme  for  solving  the 
Navier-Stokes  equations  is  fairly  simple,  employing  two-point  upwind  differ¬ 
encing  for  the  advective  terms.  Nevertheless,  the  VORTEX  code  generated 
pressure  information  that  would  have  been  quite  useful  in  trouble-shooting  and 
predicting  the  onset  of  cavitation  in  the  Taylorsville  outlet  works.  Since 
there  is  nothing  unusual  about  the  Taylorsville  configuration,  it  is  likely 
that  the  code  will  do  about  as  well  for  other  structures.  Further  study  and 
application  are  needed  to  establish  error  bounds  and  code  utility  for 
hydraulic  structures  in  general. 
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Figure  *4.  Physical-model  section  and  boundary-fitted  finite 
difference  grid  for  flow  simulation  through  flood-control 
facility  with  both  gates  open  100  percent 
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Figure  5.  Physical-model  section  and  boundary-fitted 
finite  difference  grid  for  flow  simulation  through 
water-quality  facility  with  emergency  gate  closed 


Gate  Opening  =  100  percent 
Pool  Elevation  =  570.5  ft 


Flow  Rate  =  5,000  cfs 


Figure  8.  Streamlines  for  flow  through  flood-control  facility  with 

both  gates  open  100  percent 
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Figure  9.  Comparison  of  computed  and  measured  pressure  distributions 
for  flood-control  flow,  both  gates  open  100  percent,  nodes  1-50 
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Figure  11.  Comparison  of  computed  and  measured  pressure  distributions 
for  flood-control  flow,  both  gates  open  100  percent,  nodes  4-20 


Cate  Opening  =  25 
Pool  Elevation  =  59 
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Figure  12  Velocity  vectors  for  flow  through  water  quality  facility 
Figure  12.  wiyth  service  gate  open  25  percent 
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Figure  14.  Comparison  of  computed  and  measured  pressure  distributions 
for  water  quality  flow  with  service  gate  open  25  percent,  nodes  10-40 
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Figure  15.  Comparison  of  computed  and  measured  pressure  distributions 
for  water  quality  flow  with  service  gate  open  25  percent,  nodes  2-12 
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Figure  16.  Comparison  of  computed  and  measured  pressure  distributions 
for  water  quality  flow  with  service  gate  open  25  percent,  nodes  4-16 


Gate  Opening  =  50  percent 
Pool  Elevation  =  530  ft 


Flow  Rate  -  1 ,520  cfs 


VYWV  vuwvm 

V\\\\\\\\\\\W 
"'W\\\\\\\\ 


s\v\"« 


. 

.  ,...««»* 

- NNN‘\\\\\\\\\\\W^n'"" .  ,.t|i 


*XNW\\\\\\ 
'N'W\\\\\X 
,SNVAW\\\\' 
;^w\w\vN 


.III' 

.,>1 


\\ 

XWVs'^'"'"" 


Gate  Opening  =  50  percent 
Pool  Elevation  =  530  ft 


Figure  18.  Streamlines  for  flow  through  water  quality  facility 
with  service  gate  open  50  percent 
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APPENDIX  A:  THEORY  USED  IN  THE  VORTEX  CODE 

1 .  For  Incompressible  flow,  the  governing  equations  are  the  Navier- 
Stoxes  equations  for  the  conservation  of  momentum  and  mass,  respectively: 


p£ut  ♦  V  •  (uu)]  -  yV2u  -  Vp  +  pg 


V  •  u  »  0  (A2) 

In  two  dimensions,  the  continuity  equation  (A2)  is  readily  satisfied  when  the 
components  of  u  are  derivatives  of  the  stream  function  V  : 


u  -  H 


v  =  -¥ 


The  pressure  can  be  eliminated  from  the  governing  equations  by  taking  the  curl 
of  the  momentum  equation  (A1)  and  replacing  the  two-dimensional  (2-D)  Navier- 
Stokes  equations  by  the  new  system 

p[ct  ♦  (uc)x  +  (vc)yj  -  i iV2c  ♦  (p82)  ~  (p^)  (A5) 

x  y 

v2f  -  -c  (A6) 

Equation  A5  replaces  Equation  A1 ,  subject  to  the  Boussinesq  approximation*; 
and  Equation  A6  expresses  the  relation  between  the  vorticity  and  the  deriva¬ 
tives  of  the  velocity  components 

C  -  vx  -  uy  (A7) 

subject  to  Equations  A3  and  A*4 .  The  new  system  of  equations  represents  the 
Navier-Stokes  equations  in  stream-functlon/vorticity  form,  which  is  often  the 
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most  convenient  and  manageable  form  for  calculating  2-D  internal  flow.  After 
solving  these  equations  for  V  and  x,  ,  the  velocity  components  can  be  found 
from  Equations  A3  and  AM;  and  the  pressure  can  be  calculated  by  integrating 
the  x-  and  y-components  of  the  momentum  equation: 


p[ut  +  (u2)x  +  (UV)yJ  =  HP2,J  "  Px  + 
p[vt  +  (uv)x  +  (v2)y]  =  uV2v  -  Py  +  pg2 


Equations  A8  and  A9  can  be  put  in  more  tractable  form  by  using  the  2-D  con¬ 
tinuity  equation 


ux  +  vy  =  ° 


(A10) 


along  with  Equation  A7  to  obtain 


Px  =  ~P(ut  +  uux  +  vvx)  +  PVG  _  PCy  +  P81 


Py  =  -P(vt  +  uuy  +  Wy)  -  put  +  ucx  +  pg2 


(All  ) 


(A12) 


Equations  All  and  A12  now  combine  to  form  a  single  equation  for  the  pressure 


increment 


dp  1  2  2.  ,  .  . 

-^  =  ^d(u  +  v  )  +  (g,  +  vc  -  vt  -  u.  )  dx 
P  c  I  y  t 


+  (g2  -  u c  ♦  vcx  -  vt)  dy 


(A13) 


This  final  form  is  particularly  convenient  for  Integration,  because  dp  is  an 
exact  differential.  As  a  result,  the  pressure  change  can  be  calculated  along 
any  path  defined  by  incremental  displacements  dx  and  dy  . 

2.  The  stream-function/vorticity  equations  are  solved  numerically  by 
the  VORTEX  computer  code,  which  has  been  developed  from  the  WESSEL  code 
(Thompson  and  Bernard  1985*).  WESSEL  solves  Equations  A8,  A9,  and  A10, 
whereas  VORTEX  solves  Equations  A5  and  A6.  As  a  rule,  the  latter  equations 

*  See  References  at  end  of  main  text. 


LNMivlVLv: 


are  easier  and  less  expensive  to  solve  when  using  finite  differences.  The 
input/output  structure  of  VORTEX  is  nearly  identical  with  that  of  WESSEL,  and 
both  codes  use  boundary-fitted  finite  difference  grids  generated  by  the  WESCOR 
computer  code  (Thompson  1983).  Additional  documentation  for  the  VORTEX  code 
will  be  provided  in  the  near  future. 


APPENDIX  B:  NOTATION 


Incremental  operator 

Gravitational  acceleration 

Gravity  vector 

x-  and  y- components  of  g  , 

Pressure 

Time 

Velocity  vector 

x-  and  y-components  of  u  , 

Cartesian  coordinates 

Vortici ty 

Dynamic  viscosity 

Kinematic  viscosity 

Density 

Stream  function 
Gradient  operator 
Laplacian  operator 


respectively 


respectively 
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